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ABSTRACT 



o 



, Context. Low mass starless cores present an inhomogeneous chemical composition. Species like CO and CS deplete at their dense 
I/") ' interiors, while N2H + and NH3 survive in the gas phase. As molecular line observations are used to determine the physical 
' conditions and kinematics of the core gas, chemical inhomogeneities can introduce a serious bias. 
Aims. We have carried out a molecular survey towards two starless cores, L1498 and L1517B. These cores have been selected 
for their relative isolation and close-to-round shape. They have been observed in a number of lines of 13 molecular species in 
pH , order to determine a self-consistent set of abundance profiles. 

O I Methods. In a previous paper we modeled the physical structure of L1498 and L1517B. Here we use this work together with a 
spherically-symmetric Monte Carlo radiative transfer code to determine the radial profile of abundance for each species in the 
survey. Our model aims to fit simultaneously the radial profile of integrated intensity and the emerging spectrum from the core 
center. 

^> Results. L1498 and L1517B present similar abundance patterns, with most species suffering a significant drop toward the core 
• • center. This occurs for CO, CS, CH3OH, SO, C 3 H 2 , HC 3 N, C 2 S, HCN, H 2 CO, HCO+, and DCO+, which we fit with profiles 
t ^ having a sharp central hole. The size of this hole varies with molecule: DCO + , HCN, and HC3N have the smallest holes, while 
, SO, C2S and CO have the largest holes. Only N 2 H + and NH3 seem present in the gas phase at the core centers. 
j_j ■ Conclusions. From the different behavior of molecules, we select SO, C2S, and CH3OH as the most sensitive tracers of molec- 
' ular depletion. Comparing our abundance determinations with the predictions from current chemical models we find order of 
magnitude discrepancies. Finally, we show how the "contribution function" can be used to study the formation of line profiles 
from the different regions of a core. 

Key words. ISM: abundances - ISM: clouds -ISM: molecules - stars: formation - ISM: individual(L1498, L1517B) 



1. Introduction 

Recent observations and modelling of low-mass star- 
less cores show that a number of molecular species 
deplete from the gas phase at the den se interior of 
these simplest star- fo rming r egions llKuiper et al. I 
1996t IWiilacv et al. I ll99St iKramer et al. I Jl99S : 
Alves et al. lll999UCa.selli et al. lll999HBergin et al. iboOlt 
Tafalla et al. Il2002t iBacmann et al. Il2002t IPaeani et al. 1 
20051) . The depletion of molecules in the pre-stellar 
material has important consequences for the study of the 
initial conditions of stellar birth. Molecular emission is 
routinely used to trace the physical properties and the 
kinematics of the star-forming gas, so the removal of 
certain species from the gas phase introduces a potential 
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distortion in molecular line observations. Depletion, in 
addition, systematically increases with time and gas 
density. Thus, if properly understood, it can provide a 
reliable clock to time the contraction history of dense 
cores. Understanding the chemical composition of starless 
cores has therefore become necessary to reconstruct the 
earliest phases of low-mass star formation. 

Previous studies of core chemical composition show 
that molecular depletion is a highly selective process. 
Molecules like CO and CS disappear rapidly from the gas 
phase, while species like N2H + and NH3 survive much 
longer at high densities (see references above) . As a result 
of this behavior, a core gradually develops a differenti- 
ated interior characterized by a center rich in depletion- 
resistant species surrounded by layers richer in depletion- 
sensitive molecules. Understanding this abundance pat- 
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tern is critical to interpret molecular line observations of 
cores because different species will systematically trace 
different layers depending both on their response to de- 
pletion and their level excitation. 

Chemical models already provide an indication of how 
the different molecular species will be distributed inside a 
core, but systematic molecular surveys are still needed to 
obtain a realistic picture of the core chemical composition. 
Such surveys should be carried out toward starless cores of 
simple geometry, so the abundances of the different species 
can be unambiguously derived from observations. In this 
paper, we present a survey toward two Taurus-Auriga 
cores, LI 498 and L1517B. These two cores are reasonably 
isolated and present close to round shapes when observed 
in the millimeter continuum, so they constitute id eal tar- 
gets for a systematic study. In lTafalla et al. I l|2004(l (paper 
1 hereafter), we used millimeter continuum data together 
with C 18 0, CS, N 2 H+, and NH 3 line observations of these 
cores to derive their radial profiles of density, tempera- 
ture, turbulent linewidth, and line-of-sight velocity. Now 
we complement this study with additional observations of 
a number of molecules, and use the already derived phys- 
ical models of the cores to determine in a self consistent 
manner the radial abundance distribution of 13 molecular 
species. 

In the following sections, we present the details of 
our radiative transfer modeling of the observed lines and 
the set of abundance profiles derived from this analysis 
(section 4). Using these abundances, we study the differ- 
ences and similarities between the two cores, as well as 
their chemical relation with other cores of well determined 
abundances (sections 5.1-5.3). We also use our data to test 
core chemical mode ls, in particular the recent work by 
lAikawa et al. I l|2005jl (section 5.4), and to study the differ- 
ent sensitivity of molecules to depletion (section 5.5). We 
finish using our chemical determinations and the radiative 
transfer modeling to study how the different molecular 
lines originate at the core interior and therefore reflect (or 
miss) its internal structure when used to trace the core 
gas (section 5.6). For this analysis we use as a tool the 
contribution function. 

2. Observations 

We observed L1498 and L1517B with the IRAM 30m tele- 
scope during several runs between 1999 October and 2002 
November. We made maps of these cores in the lines shown 
in Table 1, always observing in frequency switching mode 
with several receivers simultaneously and a typical sam- 
pling of 20". As a backend, we used a facility autocorrela- 
tor that provided a typical velocity resolution of 0.03-0.04 
km s _1 . During the observations, the telescope pointing 
was corrected making frequent cross scans of bright con- 
tinuum sources, and the atmospheric attenuation was cal- 
ibrated observing ambient and liquid nitrogen loads. The 
telescope intensity scale was converted into main beam 
brightness temperature using standard efficiencies. The 



Table 1. Observed lines and rest frequencies used in this 
work. 



Line 


Telescope 


Frequency 


Et 


Ref. 






(MHz) 


(K) 




CH 3 OH(J fc =2 -lo)A + 


IRAM 


96741.375 


2.3 


(1) 


CH 3 OH(J fc =3 -2 )A+ 


IRAM 


145103.200 


6.7 


(2) 


CH 3 OH(J fe =2-i-l-i)E 


IRAM 


96739.362 


0.0 


(1) 


CH 3 OH(J fe =3_i-2_i)E 


IRAM 


145097.450 


4.6 


(2) 


SO(JN=32-21) 


IRAM 


99299.890 


4.5 


(2) 


SO(JN=43-32) 


IRAM 


138178.670 


9.2 


(2) 


C-C 3 H2(Jif_ 1 X +1 =2l2-loi) 


IRAM 


85338.894 


0.0 


(3) 


C-C 3 H2(Jif_ 1 if + i=3l2-22l) 


IRAM 


145089.610 


6.7 


(3) 


Hov.AJf.Jfi' l<~ — 2i1 — 111 ) 
1 L Z ' y <J 1\ 1 ■"■-(- 1 iJ i *i - L ll/ 


IRAM 


140839.502 


0.0 


(4) 


H2CO(Jx_ 1 k +1 =2n-lio) 


IRAM 


150498.334 


0.3 


(5) 


HC 3 N(J,F=4,5-3,4) 


100m 


36392.363 


2.6 


(4) 


HC 3 N(J,F=10,11-9,10) 


IRAM 


90979.002 


19.7 


(4) 


HCO+(J=1-0) 


FCRAO 


89188.523 


0.0 


(5) 


HCO+(J=3-2) 


IRAM 


267557.619 


12.9 


(5) 


H 13 CO + (J=1-0) 


FCRAO 


86754.288 


0.0 


(6) 


DCO + (J=3-2) 


IRAM 


216112.582 


10.4 


(7) 


C 2 S(JN=67-56) 


IRAM 


86181.391 


19.2 


(8) 


HCN(JF=12-01) 


FCRAO 


88631.847 


0.0 


(4) 


HCN(JF=34-23) 


IRAM 


265886.487 


12.8 


(4) 


H 13 CN(JF=12-01) 


IRAM 


86340.168 


0.0 


(4) 



References: (1) IMiiller et all J2004I) : (2) this work; (3) 
C. Gottlieb, priv. c omm. ; (4) CDMS; (5) JPL; (6 ) 
ISchmid-Burgk et al. I |2004|) : (7) ICaselli fc Dorel fc005ft : 
f8) lYamamoto et al. I il990h 



FWHM of the telescope beam varied with frequency from 
27.7" at 90 GHz to 11" at 230 GHz. 

We observed L1498 and L1517B in HCO+(1~0), 
H 13 CO+(1-0), and HCN(l-O) with the FCRAO 14m tele- 
scope 1 in 2001 April. We used the QUARRY array receiver 
in frequency switching mode together with the facility cor- 
relator, which provided a velocity resolution between 0.03 
and 0.07 km s _1 . Observations of SiO masers were used 
to correct the telescope pointing, and an ambient load 
was used to calibrate the atmospheric attenuation. An ef- 
ficiency of 0.55 was used to convert the telescope units into 
mean beam brightness temperature. The typical sampling 
of the maps was 30", and the FWHM of the telescope 
beam at 86-89 GHz was approximately 55". 

For the detailed line modeling presented here, ac- 
curate rest frequencies are required. We have searched 
for such frequencies using on-line compilations like the 
Cologne Database for Molecular Spectros copy (CDMS, 
I Miiller et al"1l200lj) and the JPL Catalog {Pickett et al. I 
199&). In most cases, the CDMS and JPL frequencies are 
consistent with each other, so we have chosen the one 
quoted as having smaller uncertainty (where possible, we 

1 FCRAO is supported in part by the National Science 
Foundation under grant AST 94-20159, and is operated 
with permission of the Metropolitan District Commission, 
Commonwealth of Massachusetts. 
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have referred to the original determination). For HCO + , 
however, the values in the two catalogs are inconsistent, 
probably due to the difficulty measuring frequencies in 
ions, and we have preferred the JPL value for its bet- 
ter agreement with our line data. For some transitions 
of CH3OH and SO, no accurate frequencies were found, 
and we have determined them by fitting our L1498 spec- 
tra assuming an LSR velocity of 7.80 km s _1 (as mea- 
sured from the other lines). Fortunately, these few spectra 
present gaussian lines, so the astronomical determination 
is likely to be accurate (about 20 kHz). A summary of the 
frequencies used in this work is presented in Table 1. 

3. Overview of the molecular emission 

In Figures 1 and 2 we present the integrated intensity 
maps of all the lines observed toward L1498 and L1517B 
together with maps of the dust continuum emission and 
N2H + (l-0) already analyzed in paper 1. For each line, 
the integrated intensity map reflects the combined effect 
of excitation, optical depth, and molecular abundance, so 
its interpretation requires detailed radiative transfer mod- 
eling. Even without such an analysis we can appreciate the 
need for strong abundance variations by noting that most 
lines are not very optically thick and that their excitation 
increases with density toward the dust continuum peak. 
Thus, if the abundance of a molecule were spatially con- 
stant, its emission map would present a well-defined max- 
imum at the dust peak. Although this occurs for N2H+ 
(whose abundance is close to constant, paper 1), it is not 
the case for the rest of the species. 

In the larger and better resolved L1498 core, the maps 
of C 3 H 2 , H 2 CO, CH3OH, SO, HCO+ (plus isotopes), HCN 
(plus isotope), HC3N, and C2S all present distributions 
that differ from the centrally concentrated dust or N2H" 1 ". 
In a few species, like CH 3 OH, SO, and H 2 CO, the emission 
forms an almost-complete ring around the dust peak. In 
others, the emission presents discrete peaks, especially to 
SW and a NE of the dust peak, reminiscent of a broken 
ring. Such ring-like distributions are also seen in C 18 
and CS, and in paper 1 it was shown that t hey reflect the 
presence of a central deple tion hole (see also lKuiper et al. 1 
199(4 IWillacv et al. lll998|) . The maps in Fig. 1 show now 
that depletion holes must be the rule for most species in 
L1498. 

For the more compact L1517B core, the pattern of line 
emission is similar to that in L1498, although less distinct 
for some species because of angular resolution. In this core, 
most molecules present a single emission peak offset to the 
west from the dust/N2H + peak, being the most extreme 
example that of SO. This pattern is again similar to that 
found in C 18 and CS (paper 1), and for the same reasons 
as in L1498, it requires a central depletion hole. 

A simple inspection of the maps in Figs. 1 and 2 shows 
that different molecular species have depletion holes of 
different sizes. Using again the larger L1498 core as a ref- 
erence, we note that CH3OH presents a rather prominent 
hole, while DCO + has a central depression barely resolved 



by our observations. This variety of hole sizes most likely 
arises from the different sensitivity to depletion of the dif- 
ferent species and from the complex chemical changes re- 
sulting from the freeze out of CO, the main source of car- 
bon in the gas phase. The freeze out of CO decreases the 
abundance of the different C-bearing species by different 
degree, while it enhances the abundanc e of deuterium- 
bearing species ((Dalgarno fc Lenr)"lll984h . In Section 5.2 
we will present a quantitative comparison of hole sizes us- 
ing the results of a radiative transfer analysis. 

Given the large effect of depletion in the maps of Figs. 

1 and 2, it seems clear that in most cases the shape of 
the emission reflects more the chemical composition of a 
core than its true physical structure. As chemical inhomo- 
geneities seem the rule in starless cores, this high sensitiv- 
ity to chemistry of the maps should be carefully consid- 
ered when deriving physical properties from observations 
of line emission. The risk of overlooking it is graphically 
illustrated by the maps of L1498 in species like CH3OH, 
SO, or C3H2. These maps show highly fragmented dis- 
tributions with multiple peaks along an approximately 
elliptical shell, and a naive interpretation of the emis- 
sion peaks as distinct physical structures would lead to 
a picture of a highly clumped core. This is of course in 
contradiction with the the distribution of the most reli- 
able tracers (1.2mm dust continuum, N2H + , and NH3), 
which shows that the core is smooth and centrally con- 
centrated. Any correspondence between map peaks and 
core substructure, therefore, requires careful consideration 
core chemistry and self consistency check using multiple 
species. The traditional warning against the use of op- 
tically thick tracers to infer physical properties of cores 
should now be expanded to avoid using depletion-sensitive 
species for the same purpose. Comparing in Figs. 1 and 

2 the similar maps of thick and thin tracers like HCO + 
and H 13 CO+ (also HCN and H 13 CN) with the distinct 
distributions of N2H + and the dust continuum, we see 
that the danger of using depletion-sensitive molecules can 
sometimes exceed that of using thick or even self-absorbed 
tracers. 

4. Abundance analysis 

To improve the qualitative abundance analysis of the pre- 
vious section we need to model the emission of the ob- 
served species. Modeling this emission requires first de- 
termining the physical structure of each core and then 
following the transfer of radiation. For the first step, we 
make a physical model of the core that describes its den- 
sity, temperature, and gas motions, and we do so assuming 
spherical symmetry because of the close-to-round shape of 
the continuum maps. Once the core is modeled, it becomes 
like a laboratory of known physical conditions where the 
line intensities can be inverted into abundance estimates. 
For this second step, we use a non-LTE Monte Carlo code 
that solves numerically the transfer of the line emission. 
The details of these two steps are described in the next 
section. 
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Fig. 1. Integrated intensity maps from our molecular survey of L1498. The two top left panels show the 1.2 millimeter 
continuum and N2H + (l-0) data presented in paper 1. In contrast with these centrally concentrated distributions, the 
other molecular-line maps are ring-like or centrally-depressed, and have relative minima at the 1.2mm/N2H + peak. 
This is an indication of a central drop in the abundance of most species (see text). Both IRAM 30m and FCRAO 
maps have been convolved with a gaussian beam of 20-30 arcsec to filter out high-frequency noise, and for all lines, 
the contours are at 15, 30, ... percent of the peak intensity (peak intensity in K km s _1 is indicated on the top right 
corner of each panel). To better compare the survey line maps with the N2H + (l-0) emission, the contour at 75% of the 
peak is superimposed in red. The dots indicate the original map sampling. Central coordinates are aigso = 4 h 7 m 50?0, 
<5 1950 = 25°02T3'.'0. 



4.1. Core physical model and Monte Carlo radiation 
transfer 

The physical models of LI 498 and L1517B were derived 
in paper 1 from the simultaneous fit of the dust contin- 
uum, C 18 0/C 17 0, CS/C 34 S, N 2 H+, and NH 3 emissions. 
These data constrain the core density, temperature, and 



kinematics, and here we use the same parameterization for 
consistency. As mentioned in paper 1, we search for the 
simplest expressions consistent with the data. We select 
the core center from the continuum emission, and fit the 
radial profile of this emission assuming a density profile of 
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Fig. 2. Same as Fig 1 but for the L1517B core. Central coordinates are ai 950 = 4 h 52 m 7?2, S 195Q = 30°33'18'.'0. 



the form 



n(r) 



n 



1 + {r/roY 



where no, ro, and a are the central density, half- maximum 
radius, and asymptotic power law. In this way, we esti- 
mate central densities of 0.94 10 5 and 2.2 10 5 cm -3 for 
L1498 and L1517B, respectively. (Recent density deter- 
minations for these two cores bv iKirk et al. ll2005l using 
SCUBA data seem to confirm our estimates, although 
IShirlev et al. 1 l2005l have derived also from SCUBA data a 
central density for L1498 that is more than 3 times lower 
than ours.) As the cores are embedded in clouds, we trun- 
cate the analytic density profile at r = 4 10 17 cm (190" at 
140 pc) and continue it with a constant density envelope 
of 10 3 cm -3 for another r = 4 10 17 cm. This extension 



only affects the modelling of self absorbed lines and there- 
fore has no influence on our abundance profiles, which we 
model with thin species. The temperature profile of the 
cores is described with a constant value of 10 K for L1498 
and 9.5 K for L1517B (from NH3 emission data), and the 
non-thermal linewidth is also taken as constant in the in- 
ner core, with a FWHM of 0.125 km s _1 for both cores (a 
larger value is assumed in the extended cloud, see paper 
1). Finally, the two cores are assumed to be static close to 
the center (r < 1.75 10 17 cm for L1498 and r < 1.5 10 17 
cm for L1517B) and have a slow gradient outside (inward 
for L1498 and outward for L1517B). These outer motions 
are subsonic and most likely represent the kinematics of 
the extended cloud (e.g., Appendix B). 
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To solve the radiative transfer inside each core we also 
assume spherical symmetry. In paper 1 we saw that the 
maps of C 18 and CS are not circularly symmetric, in 
contrast with the maps of more reliable tracers like the 
mm continuum, N2H + , and NH3. Figs. 1 and 2 show 
now a similar situation for the species of our survey: in 
L1498 the emission is systematically brighter towards the 
SE (some species present a secondary maximum toward 
the NW), and in L1517B the western half of the core is 
brighter than the eastern half. These distributions suggest 
that the abundance of most species is not spherically sym- 
metric, despite the symmetric gas distribution inferred for 
the two cores. As studied in paper 1, the deviations seem 
correlated with the velocity of the gas, and this can be un- 
derstood as the result of differential depletion caused by 
a non spherical contraction of the cores. 2 Modeling these 
asymmetric distributions requires a 2D or even 3D radia- 
tive transfer code, which exceeds the scope of this paper. 
In the following discussion we assume spherical symmetry 
and fit for each species a circular average of the emission. 
In this way, our abundance estimates represent azimuthal 
averages over the core and therefore correspond to the re- 
sult of a mean contraction. 

As in paper 1, we solve the radiative transfer in- 
side each core with the non-LTE Monte Carlo code of 
iBernes I <|l979|) . which we have modified to include addi- 
tional molecules (see Appendix A for a summary of the 
molecular parameters used in this work) . For each molec- 
ular species, we run the code together with the core phys- 
ical model to produce a set of expected intensity distri- 
butions for different abundance profiles. These intensity 
predictions are convolved with the appropriate gaussian 
beam and compared with the observed radial profile of in- 
tegrated intensity and the central emerging spectrum for 
as many transitions of the species as we have observed 
(usually two); the best fit model is the one that fits all 
these constraints simultaneously. As the only free param- 
eter in this process is the abundance profile, observations 
of one transition are in principle enough to constrain the 
solution. Fitting at least two transitions simultaneously, 
as we do, checks the self-consistency of the process. 

4.2. Fitting procedure 

As when fitting the core physical parameters, we aim for 
the simplest abundance profiles consistent with the data. 
For each species, we start with a constant abundance 
model that fits the integrated intensity at a fiducial outer 
radius. This radius is chosen as 75" for L1498 and 55" for 
L1517B and represents a compromise between the need 
of a large enough radius (to detect central molecular de- 
pletion) and the need of a bright enough signal (to make 
a reliable fit). As shown by the dashed lines in Figures 
3 and 4, the constant abundance models systematically 
overestimate the intensity toward the core center. 

2 Note that the caption of Fig. 13 in paper 1 erroneously 
labels the N2H 4 " emission in L1498 as red when is blue. 



To improve the fit, we decrease the abundance toward 
the core center by introducing a step- function at rv The 
abundance inside r% is taken to be negligible (10~ 4 of the 
outer abundance), and the value of Th is used as a free 
parameter to improve the quality of the fit. This is the 
same approach used in paper 1 to derive abundance pro- 
files for C 18 and CS, and for most species, it provides a 
reasonable fit. In a few cases, a slightly different param- 
eterization is needed, and these exceptions are detailed 
below. 

Parameterizing the abundance profiles with a step 
function is clearly a simplification of the expected de- 
pletion be havior of molecules under realistic core cond i- 
tions (e.g., iBergin fc Lange7lll997t lAikawa et al. Il2005|) . 
However, given the fast abundance drops predicted by 
these models and the finite spatial resolution of our obser- 
vations, this parameterization is almost as good a descrip- 
tion, and it has the advantage of allowing an easy com- 
parison between molecules by comparing the outer abun- 
dance and rh values of their fits. In section 5 we will use 
this approach to study different aspects of the behavior of 
molecules under high density conditions. In the following 
paragraphs, we present further details on the fitting pro- 
cedure of each individual species. Readers not interested 
in such a detailed view can simply inspect Figures 3 and 
4 and move to section 5. 

CH3OH. Two transitions were observed for each of the 
A and E forms of this molecule, so we have fitted the 
two forms independently and determined E/A ratios of 1 
for both L1498 and L1517B. As Figs. 3 and 4 show, the 
constant CH3OH abundance models fail to fit the central 
intensity by about a factor of 2, while the models with a 
central hole provide reasonable fits to all data simultane- 
ously. 

SO. SO is the only non C-bearing molecule observed in 
this survey (N2H + and NH3 were studied in paper 1), so 
its behavior provides information complementary to that 
of all other species. Three SO transitions were observed to- 
ward L1498 and L1517B, but the SO(JN=12-21) line was 
not detected in either core. Our fits therefore only take 
into account the emission from SO(32-21) and SO(43-32), 
as the 12-21 non detection is trivially fitted by all mod- 
els. Again, the constant abundance models fail to fit the 
observed SO intensity toward the core centers by a wide 
margin. Models with a central hole, on the other hand, 
reproduce the intensities both at large and small radius, 
and show that this species is one of the most sensitive 
indicators of molecular depletion (section 5.5). 

HC 3 N. In L1498, the model with constant HC 3 N abun- 
dance does not reproduce the observed combination of 
compact emission and a central hole. Such a pattern re- 
quires that the HC3N abundance increases inwards at in- 
termediate radii, and that it has a sharp depletion hole 
near the core center. We fit this behavior by introduc- 
ing an inward jump in the abundance by a factor of 10 
at r = 1.8 10 17 cm (85" at 140 pc) together with the 
step central hole used for the other species. This type 
of profile fits well both the J=4-3 and 10-9 data (Fig. 
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Fig. 3. Radial profiles of integrated intensity and central spectra for all line data observed toward L1498. The squares 
in the radial profiles and the histograms in the spectra are observed data while the lines represent Monte Carlo model 
results. The blue dashed lines are models that assume a constant abundance chosen to fit the data at a relatively large 
fiducial radius (see text), and the red solid lines are the prediction from the best fit model. Note how the constant 
abundance models over predict the central intensity typically by a factor of 2 or more. The integrated intensities in 
the radial profiles (left panels) have a typical 1-a uncertainty of 0.01-0.02 K km s _1 for the IRAM data (with the 
exception of the high frequency lines HCN(3-2) and HCO + (3-2), which have 1-a uncertainties of 0.06 and 0.09 K km 
s" 1 ), 0.02-0.065 K km s" 1 for the FCRAO data, and 0.09 K km s" 1 for the 100m data. These values are smaller than 
the internal point-to-point scatter. The spectra in the right panels have been generated from the average of typically 
5 spectra within a 20" radius from the core center (50" for FCRAO data), in order to produce a relatively high S/N 
observation. To model this observation, we have generated synthetic spectra for each of the (typically 5) observed 
radial locations, and we have averaged them as we have done with the data. 
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Fig. 4. Same as Fig. 3 but for the L1517B core. 



3), together with the central 12-11 spectrum (not shown). 
It al so agrees w i th the prediction of the chemical model 
bv lRuffle et al. I l|l997|) . who find a late-time HC 3 N peak 
caused by CO depletion. In L1517B the situation is less 
clear due to the smaller size of this core. A constant abun- 
dance model barely fits the data and is inconsistent with 
the lopsided emission seen in the 4-3 map. To improve 



the fit, we have introduced a central hole with a small 
radius (r = 4 10 16 cm = 20") and a factor of 2 outer 
abundance drop at r = 1.25 10 17 cm (60"). This small in- 
ner hole, at the limit of our resolution, is consistent with 
the smaller-than-average hole found in L1498. As will be 
further discussed in section 5.2, HC3N seems to survive in 
the gas phase to higher densities than other species. 
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CsH 2 - In L1498, C3H2 presents the same combina- 
tion of compact emission and central hole seen in HC3N. 
A constant abundance model that reproduces the emis- 
sion at intermediate radii fails to fit the observations both 
at large and small radii (Fig. 3). As with HC3N, we in- 
troduce a factor of 10 inward jump in the abundance at 
r = 1.8 10 17 cm together with a central hole. This model 
fits well the emission, including the slightly self-absorbed 
2i2-loi spectrum towards the core center (Fig. 3). We 
note that although there are no the oretical predictions for 
a late-time enhancement of C3H2 l|Rufne et al. Ill99^ do 
not present results for this sp ecies), the chemis try of C3H2 
and HC3N seem related fe.g. JCox et al. Il989j) . It is there- 
fore possible that a single process explains the observed 
behavior of the C3H2 and HC3N in L1498. Unfortunately, 
only one C3H2 transition was observed toward L1517B. 
Lacking the brighter 2i2-loi transition, we have fitted the 
C3H2 abundance in this core with only a central hole and 
no abundance jump at intermediate radii. Observations of 
additional lines of this molecule are needed to clarify the 
behavior of C3H2 in L1517B. 

C2S. Only the JN=67-56 transition was observed to- 
ward both cores and, as Figures 3 and 4 show, its ra- 
dial distribution is inconsistent with a constant abun- 
dance profile. A model with a central abundance hole fits 
the emission both at large and small radii in both cores 
and shows that, like SO, C2S is highly sensitive to de- 
pletion. Our model for L1498 also fits nicely the emission 
of the JN=21-10 (22 GHz ), 43-32 (45 GHz) and 87-76 
(94 GHz) lines observed b v lWolkovitch et al.l l)l997|) . who 
with lKuiper et al.1 l)l996j) . first found central depletion for 
this molecule. A detailed comparison with these published 
data is presented in Appendix C. 

HCN and H 13 CN. The HCN(l-O) spectra show evi- 
dence for self absorption to wards both L1498 and L1517B 
(see also lSohn et al. Il2004j) . The HCN emission is there- 
fore dominated by the core outer layers and it cannot be 
used to derive an abundance profile inside the core. We 
thus base our abundance determination on the thinner 
H 13 CN(l-0) emission, which clearly shows a need for a 
central abundance drop in both cores. As in paper 1, we 
assume a 12 C/ 13 C ratio of 60, and we model the main iso- 
tope emission by scaling the H 13 CN abundance profile by 
that factor. This produces a fit of the HCN(3-2) emission 
in both cores but underestimates slightly the depth of the 
1-0 self absorptions. To improve the fit, we introduce in 
the constant density envelope that surrounds each core 
(r > 4 10 17 cm) an abundance enhancement of a factor of 
4 in L1498 and a factor of 2 in L1517B (the enhancement 
has no effect on the H 13 CN emission or even on the 3-2 
line). It should be noted that r > 4 10 17 cm corresponds 
to radii larger than 190", which are outside the region 
where our dust-based determination of the density is reli- 
able. For that reason, the outer abundance enhancements 
should be considered a fitting convenience that corrects 
the simple core+envelope model that we have used in the 
Monte Carlo radiative transfer, and not necessarily an in- 
dication of an abundance change in the outer core. It is in 



fact remarkable that such a simple parameterization can 
simultaneously fit both the radial profiles and the cen- 
tral spectra of lines with such different optical depths as 
H 13 CN(l-0), HCN(l-O), and HCN(3-2) (Figs. 3 and 4). 

H 2 GO. Our two H 2 CO lines, 2i 2 -ln and 2n-li , show 
evidence for self-absorption and, unfortunately, no thin 
isotopologue of this species was observed in the survey. 
Despite this, constant abundance models can be easily 
ruled out from the shape of the radial profiles, and a cen- 
tral abundance drop is needed to fit simultaneously the in- 
ner and outer emission. To properly model the self absorp- 
tion, we need again an abundance increase in the envelope 
(r > 190", and factors 15 and 5 for L1498 and L1517B, re- 
specti vely) , to which the sa me caveats mentioned for HCN 
apply. lYoung et al. I l|2004l) have observed two additional 
lines towards L1498 (ln-lio and 3i2-2n) and have also 
concluded that a central abundance drop is needed to fit 
the data. As shown in Appendix C, our best fit model als o 
fits the H 2 CO(3i2-2 u ) data from lYoung et al.l l|2004h . 
while it predicts a ln-lio self-absorption that is weaker 
than observed. This ln-lio absorption originates at the 
core outer envelope (Appendix C), so the failure in the fit 
results again from our simplified parameterization of the 
core outer layers. 

HCO+, H 13 CO+. Both the J=l-0 and 3-2 transitions 
of HCO + are deeply self absorbed, so our abundance de- 
termination relies on the thinner H 13 CO + (1-0) line and 
assumes again a 12 C/ 13 C ratio of 60. As in all other 
species, a central drop in the HCO + abundance is needed 
to fit the shape of the radial profiles. To fit the HCO + 
self absorption, we again need to include an outer enve- 
lope (r > 190") abundance enhancement of a factor of 4 
in the L1498 model (still the predicted 1-0 line is brighter 
than observed mostly due to poor modelling of the red 
component) while no abundance enhancement is needed 
for L1517B. 

DCO + . A central abundance hole is also required to fit 
the DCO + data, but its radius is smaller than the HCO + 
hole radius. This difference most likely results from a cen- 
tral increase in the deuterium fraction ation caused by the 
CO depletion fe.g JCaselri et al. Il999j) . which partly com- 
pensates the DCO + freeze out at the inner core. The ef- 
fect is of course incomplete because the DCO + abundance 
end ups falling at the highest densities. Still, the presence 
of DCO + inside the CO and HCO + depletion holes sug- 
gests that a small amount of CO and HCO + survives at 
high densities. Comparing the DCO + and HCO + abun- 
dances in the outer core, we derive deuteration ratios of 
0.02 and 0.03 in the outer layers of L1498 and L1517B, 
respectively. Each of these values is a factor of 2 lower 
than th e deuteration ratios measured toward the inner 
core bv ICrapsi et al. I l|2005|) using N 2 H+ and N 2 D+, in 
good agreement with the expectation of a degree of deuter- 
ation that increases toward the core center. 

Although the abundance enhancements at large radii 
required to fit the self absorptions in HCN, H2CO, and 
HCO + need confirmation using a model of the large-scale 
structure of the cores, we note that such enhancements 
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are not totally unexp ecte d. Absorption line stu dies by 
iLucas k Lisztl l|l996|) and iLiszt k Lucasl i|200l|) reveal 
that species like H 2 CO, HCO + , and HCN present large 
column densities in diffuse clouds (but not N2H + , for ex- 
ample). It is therefore possible that the lower density gas 
surrounding LI 498 and L1517B has a chemical compo- 
sition similar to the diffuse clouds studied by Liszt & 
Lucas. This would naturally explain the large amounts 
of low-excitation gas needed to produce the observed self 
absorptions. 



4.3. Ortho-to-para and isotopic ratios: total 
abundances 

For species like NH 3 , H 2 CO, and C3H2, our observations 
and modeling only determine the abundance of the ortho 
or para form of the molecule, and we need to assume an 
ortho-to-para ratio (OPR) to convert our estimate into 
a total molecular abundance. This OPR depends on the 
formation history of the molecule, and is therefore some- 
what uncertain. Here we assume that the three species are 
formed by gas-phase reactions, because grain production 
will require a mechanism to release the products to the gas 
phase, a difficult task given the strong freeze out observed. 
(This of course does not imply that NH3, H2CO, and C3H2 
do not form on dust grains, but that the gas-phase chem- 
istry is separate from the dust-grain produ ction, as seen 
in the models of, e.g.. lAikawa et al~l 1200.4 ) If this is the 
case, the OPR should be close to the high-temperature 
limit, because the energy difference between ortho and 
para forms (2.4 K for C 3 H 2 , 15.2 K for H 2 CO, and 23.4 
K for NH3) is significantly lower than the energy release d 
during molecule formation Ce.g.. lDickens k Irvine lll99^ . 
Indeed, OPRs close to 3 (high temperature limit) have 
been found for H 2 CO in s t arless cores, includ i ng L14 98, by 
iDickens k Irvine I l)l999|) . iTakakuwa et al. I l|200l() mea- 
sure a C 3 H 2 OPR of 2.4 towards TMC1, also close to 
the high temperature limit. We thus assume OPRs of 3 
for H 2 CO and C 3 H 2 and 1 for NH 3 . 

In addition to OPRs, we need to assume isotopic ra- 
tios to convert rare isotopologue abundances into main 
species values. Along this paper and in paper 1 we have 
used the standard isotopic ratios of 12 C/ 13 C = 60, 32 S/ 34 S 
= 22, and 18 0/ 17 = 3.65, to which we now add the ter- 
restrial 16 0/ 18 ratio of 500. Using these values, we esti- 
mate the final main isotopologue abundances presented in 
Table 2. In the following discussion we will refer only to 
the main molecular forms, although it should be remem- 
bered that most abundances are based on the rare (and 
optically thin) isotopologues. 

Finally, we stress the dependence of our abundance 
estimates on the assumed dust parameters. As mentioned 
in paper 1, the values of the dust te mperature and emis- 
sivity are still poorl y known (e.g.. iBianchi et al. 1 ^003: 
iKramer et aT~ll2003j) . and their uncertainty propagates 
into the determination of the gas density profile and total 
core column density. These determinations, in turn, affect 



Table 2. Molecular abundances with respect to H 2 
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7.5 10 
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3.0 10~ 9 


1.0 


3.0 10~ 9 


1.15 




N 2 H+ 


1.7 10~ 10 





1.5 lO" 10 





(1) 


NH 3 


2.8 10~ 8 


(3 = 3 


3.4 10" 8 


(3 = 1 


(2) 


CH3OH 


6.0 lO" 10 


1.2 


6.0 10 -10 


0.8 


(3) 


SO 


4.0 10~ 10 


1.2 


2.0 10~ 10 


0.9 




C3H2 


1.6 10~ 9 


1.0 


9.3 10~ 10 


0.6 


(4) 


HC 3 N 


1.0 10~ 9 


0.85 


4.5 10~ 10 


0.4 


(5) 


H 2 CO 


1.3 10~ 9 


1.3 


6.7 10" 10 


0.7 


(6) 


HCO+ 


3.0 10~ 9 


1.15 


1.5 10~ 9 


0.6 


(7) 


DCO+ 


5.0 lO" 11 


0.65 


5.0 10" 11 


0.6 




C 2 S 


4.0 10 -10 


1.25 


1.0 10~ 10 


0.8 




HCN 


7.0 10~ 9 


0.8 


3.0 10~ 9 


0.53 


(8) 



(1) Factor of 3 drop for r > 1.8 x 10 17 cm in L1498 

(2) OPR = 1 assumed, X(NH 3 ) = X {n{r)/n f 

(3) E/A = 1 

(4) OPR = 3 assumed and factor of 10 drop for r > 1.8 x 10 17 
cm in L1498 

(5) Factors of 10 (L1498) and 2 (L1517B) drops for r > 1.8 x 
10 17 cm 

(6) OPR = 3 assumed. Factors of 15 (L1498) and 5 (L1517B) 
enhancement for r > 4 x 10 17 cm 

(7) Factor of 4 enhancement for r > 4 x 10 17 cm in L1498 

(8) Factors of 4 (L1498) and 2 (L1517B) enhancement for r > 
4 x 10 17 cm 



the abundance estimates through their combined effect on 
the level excitation and the H 2 column density. For most 
species, the abundance estimate uses optically thin lines 
from low-energy levels, so the effect of the density on the 
excitation is smaller than its effect on the column density. 
Thus, to first order, we can correct for dust parameter 
changes with the following simple equation: 



X(K,T d ) = A(0.005, 10) 



0.005 g cm" 



B v {T d ) \ 
B U (10K)J ' 



where B v is the Planck function, A(0.005, 10) is our 
abundance determination using a 1.2mm dust emissiv- 
ity Ki,2mm = 0.005 cm 2 g _1 and a dust temperature of 
Td — 10 K, and X(k,Tj) is the corresponding abundance 
for arbitrary values of re i .2mm. and T^. If, for exam ple, 
Ki.2mm = 0.01 cm 2 g" 1 jO ssenkopf k He nning Il994j) and 
T d = 8K l)Evans et al. l200lHGalli et al. 12002^1 . our abun- 



dance values will need to be multiplied by 1.4. 
5. Discussion 

5.1. Comparison between L1498 and L1517B 

L1498 and L1517B present similarities and differences in 
their physical properties. They have almost the same ki- 
netic temperature, level of turbulence, and total gas col- 
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Fig. 5. Abundance comparison between L1498 and 
L1517B. Top: ratio between the abundances outside the 
central hole (Xo) in the two cores (for HC3N and C3H2, 
Xo is the abundance between the hole and the outer cut- 
off). For most species, this ratio is close to 1 (within a 
factor of 2), indicating a similar pre- freeze out composi- 
tion. Bottom: ratio between central depletion holes (rh). 
L1498 presents systematically larger hole ratios by a fac- 
tor of about 1.5, which is consistent with the the visual 
impression from the integrated maps. 



uma density, while they differ by a factor of 2 in cen- 
tral density and half maximum density radius (paper 1). 
Analogously, the two cores present similarities and dif- 
ferences in their chemical composition. To examine them 
here, we compare the X$ and values derived for the 
different species. 

The top panel of Figure 5 shows the ratio between the 
outer molecular abundances derived in L1498 and L1517B 
for each molecule in our survey (values for C 18 0, CS, 
N 2 H+, and NH 3 are from paper 1). As can be seen, there is 
good agreement between the X parameters estimated for 
the two cores, with most differences being smaller than 
a factor of 2. The main outliers in the plot are C2S (4 
times more abundant in L1498 than in L1517B) and C 18 
(3 times more abundant in L1517B), and although their 
abundance differences could be real, it should be noted 
that their Xq determination is specially prone to error: 
C2S is depleted up to such a large radius in L1498 that 
Xo is determined using very limited information of the 
outer core, and the Xo value for C 18 can be contami- 
nated by the extended emission from the ambient cloud. 
Given these and other uncertainties, it is rather remark- 
able that most Xo values agree within less than a factor 
of two, and we take this fact as an indication that the two 
cores share very similar molecular compositions in their 
outer (undepleted) layers. 
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Fig. 6. Normalized radius of the abundance hole rh for all 
molecules in the survey. The red squares represent L1498 
data and the blue triangles are L1517B data (normaliz- 
ing radius is 1.04 x 10 17 cm for L1498 and 0.66 x 10 17 for 
L1517B). Species like CO, SO, and C2S present larger rh 
values, while HC3N, DCO + , and HCN have smaller val- 
ues and therefore must remain in the gas phase to higher 
densities. 



The bottom panel of Fig. 5 presents a comparison be- 
tween the radius of the central abundance holes in L1498 
and L1517B. In addition to some scatter, the figure reveals 
a trend for 77, to be larger in L1498 by a factor of about 
1.5. This trend agrees with the impression from the maps 
of Figs. 1 and 2 that L1498 has a larger depletion hole than 
L1517B for most or all molecular species. Unfortunately, 
the scatter and the uncertainties in the analysis make it 
impossible to determine whether the rh ratio varies with 
molecule. 

The similar outer abundances of L1498 and L1517B 
suggest that both cores have contracted from ambient 
gas with similar initial compositions. It is unclear how- 
ever whether the different rt values result from the cores 
being at different evolutionary stages or from them hav- 
ing followed different contraction paths. L1517B is more 
centrally peaked, denser at the center, and more circularly 
symmetric than L1498, so it would seem to be at a more 
advanced evolutionary stage. If this is the case, and both 
cores are following the same contraction path, the smaller 
rh of L1517B would imply that the depletion hole in a 
core shrinks as it contrasts; this could result from the con- 
traction time scale being shorter than the freeze out time 
scale during the late stages of contraction. Alternatively, 
L1517B may have contracted faster than L1498, so its cen- 
tral gas has had less time to freeze out. Further studies 
of core composition using a larger sample of systems are 
needed to clarify this important issue. 

5.2. Comparison between molecules 

The radiative transfer results of Table 2 suggest that dif- 
ferent molecules deplete at different radii. In addition to 
the extreme behavior of N2H" 1 " and NH3, which do not 
present abundance holes, several molecules have holes that 
are systematically smaller or larger than average in both 
L1498 and L1517B. Identifying these molecules can help 
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better understand the details of depletion. It can also pro- 
vide a set of tracers to selectively sample cores at different 
depths. 

To compensate for the systematically larger holes in 
L1498, we normalize the hole radii. We do this by divid- 
ing the radius of each species by the average hole radius 
in the core, which is calculated from the mean radius of 
all species with small dispersion (i.e., we exclude CO and 
CS, see below, in addition to N2H+ and NH3). This nor- 
malized hole radius is presented in Fig. 6 for all species in 
our survey. As the figure illustrates, the hole radii span a 
range of about a factor of 2. The CO and CS estimates 
present large scatter, which could result from a real dif- 
ference between the cores or from artifacts like contami- 
nation by the extended cloud in the case of CO and the 
fact that L1498 the CS analysis is based on the optically 
thin C 34 S emission while the analysis in L1517B used the 
thicker CS line (as C 34 S was too weak). All other species, 
on the other hand, present a reasonably low scatter be- 
tween their L1498 and L1517B results. This low scatter 
allows distinguishing between species with relatively large 
and small depletion holes. In the first group we find SO, 
C2S, and CH3OH, which in both cores have larger-than- 
average abundance holes. In the second group we find 
HC3N, DCO + , and HCN, whose central holes are more 
than 50% smaller than the holes in the first group. Such 
a relatively large difference in suggests a correspond- 
ing difference in the depletion behavior of these species. 
It also shows that HC 3 N, DCO + , and HCN survive in the 
gas phase up to higher densities than most other species, 
and this makes them interesting line tracers of the mid- 
dle layers of dense cores. Recent infall sear ches in star- 
less cores l|Lee et al. ll2004USohn et al. 1120041) . do use this 
property of DCO + and HCN to penetrate deep er than us- 
ing t raditional tracers like CS or HCO + (e.g.. iLee et al. I 
l200ll) . Our L1498 and L1517B analysis, however, shows 
that HCN and DCO+ are no substitute for N 2 H+ or NH 3 , 
as they end up depleting at densities of about 10 5 cm~ 3 . 



5.3. Comparison with other starless cores 

TMC1 is the core whose chemical composition has been 
studied with most detail, and its abundances are com- 
monly used as standards for comparison with other low- 
mass star-forming regions. Current estimates of the TMC1 
abundances, however, do not consider the effect of molecu- 
lar depletion, so a comparison with our L1498 and L1517B 
estimates can only be approximate. Until a more realistic 
analysi s of TMC1 is availab le, the most complete study is 
that of lPratan et al. "I lll997h . who have presented FCRAO 
observations and analysis of more than a dozen molec- 
ular sp ecies, 8 of them com mon to our L1498/L1517B 
survey. IPratap et al. ] (|l997h normalize their abundances 
to HCO + , and provide values for three positions along 
the TMC1 filament: the cyanopolyyne peak, the ammo- 
nia peak, and the SO peak. To compare with these val- 
ues, we select the L1498 and L15I7B abundances outside 



the depletion hole and normalize them to HCO + . The re- 
sulting values agree well with those of TMC1, and in no 
case the difference exceeds one order of magnitude. The 
best agreement occurs for the cyanopolyyne peak, where 
the average ratio between the L1498 and TMCI abun- 
dances for the 8 species is 1.3 ± 1 . This position also pro- 
vides the best match to the L1517B abundances, although 
the ratio is significantly larger (2.3 ± 2) because our es- 
timate of the HCO+ in L1517B is half of that in L1498. 
Given the very different methods and the uncertainties 
in the two analysis, the agreement with the TMCI esti- 
mates seems rather good. The fact that the best match 
to the abundance in the undepleted outer parts of L1498 
and L1517B occurs for the cyanopolyyne peak is also in 
good agreement with the interpretation that this TMCI 
position is the most chemically young: of the filament 
l|Suzuki et all ll 9921 Iffirahara et mI. Ill 99:1). It also sug- 
gests that the abundanc es of TMCI are not anomalously 
high l|Howe et al. Ill99fi|) . but representative of the popu- 
lation of dense cores in the Taurus cloud. 

A different set of abundance determinations in star- 
less cores has been recently provided by the Texas group 
using a technique similar to ours: density estimates using 
continuum (SCUBA) data followed by Monte Carlo mod- 
eling; of the line emission llEvans e t al. lEo Oll ILee et al. I 
l2003|) . As a result of this work. iLe^^tHiLT l)2003[) have es- 
timated the abundance of HCO+, H 13 CO+, and DCO+ 
in L1512, L1544, and L1689 by fitting step functions 
with non zero value at the core center. Comparing the 
outer abundances estimated by these authors with those 
in L1498 and L1517B, we find reasonable (factor of 2) 
agreement for HCO + in L1512 and L1544, the two Taurus 
cores, but almost one order of magnitude difference for 
H 13 CO + and DCO + in the same objects. This result is 
somewhat surprising, especially for H 13 CO + , as our abun- 
dance values ag ree with a standa rd isotopic 12 C/ 13 C ratio 
of 60, while the lLee et all l|2003|) values imply a ratio one 
order of magnitude lower. The origin of this discrepancy 
is not clear, although it most likely results from a com- 
bination of a higher dust emissivity, l ower gas tempera - 
ture, and smaller radius coverage bv ILee et all (|2003|) . 
A similar (factor of 4) discrepancy occurs between our 
H^CO abundance d etermination in L1498 and that by 
lYoung et al. I l|2004h in the same core, despite our model 
fitting the H2CO(3i2-2n) emission from these authors 
(Appendix C). This discrepancy most likely arises from 
a different cho i ce of dust emissivity, which as discussed by 
IShirlev et al. 1 1)2005|) . makes our L1498 models differ from 
that of the Texas group by a factor of 3 in H 2 column den- 
sity. Such a systematic discrepancy between models high- 
lights the urgent need for an accurate determination of 
the dust emissivity at millimeter and submillimeter wave- 
lengths. 

Better agreement between our abundance determina- 
tions and those in o ther starless cores occur for N2H+. 
ICaselli et al. I l|2002al) have estimated an N 2 H + abundance 
of (2± 1) x 10~ 10 for a sample of 25 starless cores assuming 
virial equilibrium, and this result is in excellent agreement 
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with our determination. Also. lKeto et al. 1 l|2004j) have car- 
ried out a detailed analysis of the N 2 H + emission in 3 star- 
less cores, one of them being L1517B. They derive for this 
object a central density within a factor of 2 of ours, and 
an N2H + abundance that differs from ours by less than 
10%. 

Even if the larger abundance differences between cores 
are real, it seems that in most cases they are well within 
the factor of 3-4 range. This of course does not imply 
that starless cores form an homogeneous family. The gra- 
dient across TMC1 reveals differences in abundance that 
are unlikely to be explained simply as the result of dif- 
ferential depletion, and a number of cores are known to 
have significantly lower ab undance of certain late-time 
species like NH 3 and N 2 H+ ijSuzuki et al. Ill992l) . L1521E, 
for example, has negligible depletion of C-bearing species 
and NH3 and N 2 H + abundances o ne order of magni- 
tude lower than L1498 and L1517B llSuzuki et al. Ill992t 
iHirota et a.l~1l2002t iTafalla fc Santiago l2004h . Our molec- 
ular survey of the L1521E core (in preparation) seems to 
suggest, however, that most of chemical differences be- 
tween cores can be explained as a result of differential de- 
pletion plus time evolution of the N-bearing species, and 
that the undepleted abundances in the outer layers do not 
present extreme (more than a factor few) core to core vari- 
ations. A systematic survey of a larger sample of cores is 
still needed to confirm this result. 



5.4. Comparison with the Aikawa et al. (2005) 
chemical model 

The chemical evolution of cores as they contract un- 
der gravity has been stu died by d i fferent authors. 
iBergin fc Langer ] lll997h and ICharnleTl <|l997h were first 
to model the differential depletion of C and N-bearing 
species as a result of their freeze out onto the dust grains 
together with a lower binding energy of N 2 . More com- 
plete chemical networks coupled to realistic contraction 
physics have been used to improve on this earlier work 
jAikawa et al. I|200ll iLi et al. Il2003 lAikawa et al. Il2003t 
Lee et al. Il200c L Shematovich et al. Il2003|) . Very recently. 



Aikawa et al. 1 l)2005|) (AHRC05 hereafter) have presented 



the most up-to-date chemical model of starless core con- 
traction. These authors have followed the evolution of two 
dense cores initially having Bonnor-Ebert density profiles, 
one subject to a small over-density factor (their a = 1.1 
case, where a is the gravity-to-pressure ratio), and the 
other highly over dense (a = 4.0). As intuitively expected, 
the core with the small over-density factor loses equilib- 
rium and contracts slowly (timescale 1 Myr), while the 
highly over dense core contracts in one tenth of the time. 
Due to these very different time scales, the two cores de- 
velop significantly different chemical compositions by the 
time their central densities reach values like those of L1498 
and L1517B. The a = 1.1 core has lost most molecular 
species at the center, while the a = 4.0 core still retains a 
significant fraction of depletion-sensitive molecules in the 



gas phase. From this dichotomy, AHRC05 argue that the 
a = 1.1 model approximates the evolution of quiescent, 
heavily depleted cores like L1498 and L1517B, while the 
a = 4.0 model simulates dense, but not depleted cores like 
L1521E. 

Our molecular survey of L1498 and L1517B can be 
used to test the predictions of the different chemical mod- 
els of core contraction. We choose the AHRC05 model 
because it is the most complete model available and be- 
cause it presents predictions of both the core velocity field 
and its chemical composition. To test this model, we have 
run a series of Monte Carlo radiative transfer calculations 
using the observed density profiles of L1498 and L1517B 
together with the abundance and velocity profiles pre- 
dicted by the AHRC05 models at the time when the core 
central density reaches uq — 1.5 x 10 5 cm -3 (note that 
AHRC05 quote densities of H nuclei while we use densi- 
ties of H 2 molecules). At this time, the density profiles 
of the AHRC05 models are close to those of L1498 and 
L1517B, although there are still factor-of-2 differences be- 
tween the density profiles that should be kept in mind 
when comparing with our observations. 

Concerning the velocity field, we find that the a — 4.0 
model predicts for NH3 lines much broader than observed 
(0.38 km s _1 versus 0.2 km s _1 , paper 1), and for species 
with central depletion like C 18 and CH3OH, the model 
predicts spectra with two peaks. These peaks correspond 
to the front and the back sides of the core, and are sepa- 
rated by 0.4 km s _1 because each side is moving toward 
the center with a velocity of about 0.2 km s _1 (Fig. le 
in AHRC05). Such broad, double-peaked profiles are not 
observed in LI 498 or L1517B, and they rule out a fast 
contraction model for these cores. The a ~ 1.1 model, 
on the other hand, predicts narrower spectra, similar to 
those observed. In the absence of turbulence, this model 
predicts the correct NH3 linewidth toward the center for 
both L1498 and L1517B (0.2 km s" 1 ), but for species with 
central depletion, it predicts lines that have two peaks sep- 
arated by 0.15 km s _1 , and this is not observed. Although 
these peaks can be blended into a single feature by adding 
a small amount of turbulence, the lines are already broader 
than observed and the extra turbulence further degrades 
the fit. For the narrow-lined CH3OH emission, for exam- 
ple, the AHRC05 model without turbulence predicts lines 
that are broader than measured by 20% in the case of 
L1517B and by 70% for L1498. This disagreement may be 
exaggerated by the imperfect match between the density 
profiles of our cores and the AHRC05 models, but also il- 
lustrates the fact that in the AHRC05 a = 1.1 model, den- 
sities like those of LI 498 and L1517B are reached at rel- 
atively late stages, when the core contraction has started 
to accelerate. For this reason, the model predicts that the 
L1498/L1517B phase will only last 10% of the core time 
life, or ju st 10 5 yr, which is to o short for a typical starless 
core (e.g. JLee fc Mvers Il999l) . The relatively fast contrac- 
tion of the model is therefore not supported by the obser- 
vations. Adding a magnetic field can help to slow down the 
contraction, but this may introduce additional problems 
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Fig. 7. Correction factors needed to mak e the abundances 
predicted by the lAikawa et al. I (120051) a = 1.1 model 
reproduce the observed intensities in the L1498/L1517B 
molecular survey. The blue triangles (L1517B data) and 
the red squares (L1498 data) indicate that order of mag- 
nitude corrections are needed in a number of species. 
The comparison between model and observations has been 
done at the fiducial outer radius defined in section 4.2. 
The dashed and dotted lines indicate factor of 2 and 10 
corrections, respectively. 



with the chemistry. A long period of slow contraction, for 
example 1 Myr, w ould imply C O depletion at densities as 
low as 10 cm" " 3 llLegerlll983h . which is lower than ob- 
served. Accretion of new (undepleted) material from the 
surrounding molecular cloud may help mitigate this dif- 
ficulty. An additional problem with chemical models of 
magnetized clouds is that they predict N 2 H + to trace the 
quiescent and thermally supported c ore nucleus in con 
tracting starless clouds such as L1544 llShematovich et al~ 
20031) in disagreement with observations l)Caselli et al.' 



To test the AHRC05 abundance predictions, we now 
fix the velocity field of each core to the best fit value de- 
termined in paper 1, and we use the Monte Carlo code to 
estimate predicted radial profiles of integrated intensity 
for all species observed in L1498 and L1517B and pre- 
dicted by model a = 1.1 (no = 1.5 x 10 5 cm -3 ). If the 
AHRC05 abundance profile does not fit the observations 
at the fiducial outer radius defined in section 4.2 (75" for 
L1498 and 55" for L1517B), we estimate the global factor 
by which the model abundance needs to be multiplied to 
match the data at that radius. This factor measures the 
deviation of model from the observations, and is shown 
in Fig. 7 for both the L1498 and L1517B models. As ex- 
pected from the good agreement between the L1498 and 
L1517B abundances (section 5.1), the correction factors 
derived from the two cores agree within a factor of two. 
This shows that the corrections are almost independent of 
the exact structure of the core. 



The correction factors in Fig. 7 cover more than 
three orders of magnitude, and this indicates that some 
species are seriously over or under predicted by the model. 
Considering all the uncertainties of our fitting procedure, 
we take any correction smaller than a factor of two as a 
reasonable match between model and data. This occurs 
for C3H2 (although the model misses the observed outer 
abundance drop), HCO + (although L1498 requires a fac- 
tor of 3 correction), and CO. Correction factors between 
2 and one order of magnitude are considered "clear devi- 
ations" between model and data, and include SO, DCO + , 
C2S, CS, and N 2 H + . Interestingly enough, three of these 
species are S-bearing, and all of them are under predicted 
by a factor of a few (5-7.5 for CS, and 2-4 for SO and 
C2S). This suggests that the model overestimates the de- 
pletion of S in its initial atomic conditions by a similar 
factor. Finally, correction factors larger than one order of 
magnitude indicate a "serious deviation" between model 
and data, and occur for CH 3 OH, H 2 CO, HCN, and NH 3 . 
Given its large value, the deviation of H2CO is likely real, 
although thin isotopologue observations are needed to bet- 
ter quantify the error. On the other hand, the large under- 
production of CH3OH is not surprising given the lack of a 
known gas-phase p roduction mechanism for this molecule 
l|Luca et al. 11200^1 . As for HCN, it is possible that the 
addition of photochemistry to the model may help correct 
the overproduction of this molecule. 

The need for order-of-magnitude corrections does not 
necessarily imply that the AHRC05 model fails to repro- 
duce all aspects of the chemical structure in cores like 
L1498 and L1517B. Most abundance profiles, when scaled 
appropriately, reproduce the shape of the radial distribu- 
tion of intensity, and this suggests that the model captures 
at least part of the process of molecular depletion as the 
core contracts. To study how well the AHRC05 model pre- 
dicts the radial variation in the abundance of the different 
species, we compare the shape of the observed intensity 
profiles with the shape predicted by the model after cor- 
recting the abundances by the appropriate global factor. 
We do this by defining a "concentration factor" as the 
ratio between the intensity at the core center and the in- 
tensity at a fiducial outer radius (75" for L1498 and 55" 
for L1517B, see section 4.2). Such a factor measures how 
centrally peaked the emission is, and can be easily esti- 
mated both in the model predictions and in the data. The 
ratio between the model and data concentration factors 
will equal one in a perfect match, will be larger than one 
if the model under-predicts the central abundance drop 
(so its emission is more centrally peaked) , and will be less 
than 1 the other way around. 

Figure 8 presents the ratio of concentration factors 
(model over data) in both L1498 and L1517B for all avail- 
able species. For L1498, the AHRC05 model predicts on 
average the correct size of the central hole, and the aver- 
age model-to-data ratio is approximately 1. For L1517B, 
the model over predicts the size of the central hole, as 
would be expected due to its smaller dimension (section 
5.1), and the mean ratio of concentration factors is about 
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Fig. 8. Compari son between the c oncentration ratios pre- 
dicted by the lAikawa et al. 1 l)2005|) model and the observa- 
tions of L1498 (red squares) and L1517B (blue triangles). 
A value equal to 1 indicates that the model predicts emis- 
sion with the same central concentration as observed, a 
value larger than 1 indicates that the model is more con- 
centrated than the data (small central hole), and a value 
less than 1 indicates a model prediction flatter than the 
data (too large a hole). There is an overall (factor of 2) 
agreement between model and data, although the L1517B 
data seems systematically lower than 1. The model there- 
fore overestimates the central hole seen in L1517B. 



0.75. This systematic over prediction for LI 51 7B reflects a 
global problem modeling the core, so in the following dis- 
cussion we concentrate on the L1498 results. As the figure 
shows, SO and HCN present the largest ratio of all, in- 
dicative that the model under predicts their central abun- 
dance hole by the largest factor. This occurs in SO because 
this species has a relatively large hole, while the AHRC05 
model predicts an average value. For HCN, the data show 
a relatively small central hole, but the AHRC05 model pre- 
dicts negligible depletion at the time when the core has 
the central density of L1498 and L1517B. (Note that the 
model correctly predicts a well defined, relatively smaller 
HCN hole at later times.) DCO + , on the other hand, is 
the species with smallest ratio in the figure, indicating that 
the AHRC05 model predicts a hole larger than observed. 
This is most likely due to an under prediction of the cen- 
tral deuterium enhancement, as the model predicts the 
correct depletion hole for the main isotopologue HCO + . 
Finally, the radial behavior of the two species without cen- 
tral freeze out, N2H 4 " and NH3, is well predicted by the 
model after dividing the N2H + abundance by a factor of 
3.5 and the NH 3 abundance by a factor of 20. This gen- 
eral over prediction of N-bearing species probably results 
from an underestimate of the binding energy to grains (see 
below), while the improved behavior of NH3 arises from 
the new treatment of t he N?H + dissociative re combination 
thanks to the work bv lGenpert et al. 1 l|2004 - 

A final problem affecting m ost current chemical mod els 
(including AHRC05, but also lBergin fc Langer"1ll997|) is 
that they owe most of their success in explaining the differ- 
ent depletion behavior of C-bearing and N-bearing species 
to the assumption that CO and N2 have significantly 
different binding energies to grains . Re cent laboratory 
measurements bv lQberg et al. 1 1)2005[) and lBisschop et ahl 




Fig. 9. Effect of depletion on the different lines and tran- 
sitions observed in the survey. The y-axis represents the 
ratio between the central intensities of a constant abun- 
dance model and the best fit model for each of the lines 
observed in L1498 (red squares) and L1517B (blue trian- 
gles). Large values indicate a strong effect of the central 
depletion in the emission from the core center. 



(2006), however, show that the two binding energies dif- 
fer by less than 10%, and this may seem insufficient to 
account for the different depletion behavior of CO and 
N 2 H+ (note however that lAikawa et aT~ll200ll find differ- 
ential depletion e ven when using simil ar binding energies 
for CO and N 2 ). iFlower et al. I l|2005|) have proposed an 
alternative explanation in terms of a lower sticking coeffi- 
cient for N 2 or N. While the former alternative seems rule d 
out by laboratory measurements ijBisschop et al. l|20 06). 
the latter is still a viable solution, although no laboratory 
measurements exist yet to confirm or refute the idea. This 
remaining uncertainty in our understanding of the pro- 
cess behind the differentiation of C-bearing and N-bearing 
species -the most visible feature of core chemistry- illus- 
trates how it is still premature to use current chemical 
models to assign contraction ages to cores. A new gen- 
eration of models is still needed to fulfill the promise of 
providing core studies with a reliable chemical clock. 

5.5. How to best identify depletion (or its absence) 

The simplest use of depletion as a qualitative indicator 
of core age is the classification of cores as chemically 
evolved if they show evidence for molecular depletion and 
as chemically young if they do not. Mos t dense cores, in 
fact, suffer from severe d epletion (e.g., iBacmann et al. I 
120021 iTafallaet al~ll2002l). while only a minority seem 
unaff ected by it l|Hirota et al. IbOO^ iTafalla fc Santiago! 
2004). This suggests that young cores arc rare, cither be- 
cause they are absent from clouds or because an observa- 
tional bias limits our ability to recognize them in surveys. 
In either case, it is of interest identifying the most sensi- 
tive tracer of molecular depletion to use it for systematic 
searches of young cores. Choosing such a tracer requires 
some consideration. Sensitivity to depletion depends not 
only on the relative size of the central hole, but on the 
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response of the emission to dense gas and on possible op- 
tical depth effects. In this section, we use our radiative 
transfer analysis of L1498 and L1517B to assess the effect 
of depletion on the emission of different molecules under 
realistic core conditions, in order to find the species most 
sensitive to central depletion. 



Strictly speaking, the sensitivity to depletion depends 
not on the molecule but on its transitions. These may 
greatly differ in critical density and optical depth, espe- 
cially in molecules with complex level structure. However, 
as we will see below, the results tend to agree within a 
factor of two. To measure quantitatively the sensitivity to 
depletion of a molecular transition, we compare the inten- 
sity toward the core center predicted by our best-fit Monte 
Carlo model with the prediction from a model of constant 
abundance equal to the outer value in the best fit (solid 
and dashed lines in Figs. 3 and 4). The ratio of the con- 
stant abundance result over the depletion result measures 
how much brighter a core with constant abundance would 
appear compared with the core suffering real depletion. 



Figure 9 presents the depletion sensitivity ratio for all 
transitions observed toward L1498 (squares) and L1517B 
(triangles). For most molecules, the L1498 and L1517B 
values agree within a factor of two even when multi- 
ple transitions have been observed. As the figure shows, 
a sensitivity factor of 2 (dashed line) is common for 
most species, indicating that molecular depletion typi- 
cally halves the expected intensity towards the core center. 
Several species, however, consistently present sensitivity 
ratios larger than 2. The SO lines reach the highest values 
of all, and have a mean factor of 4. This behavior is some- 
what expected from the systematically larger depletion ra- 
dius of this molecule (Section 5.2) and its moderately large 
dipole moment (1.55 D), which makes it sensitive to the 
presence or absence of SO in the inner core. Not surpris- 
ingly, the highest SO ratios correspond to the JN=43-32 
transition (138 GHz), which has an Einstein A coefficient 
larger than the other observed SO line, JN=32-21 at 99.3 
GHz. Second in sensitivity to depletion is C2S, with small 
scatter and also a mean value of 4. This molecule again 
combines a larger-than-average depletion radius with lines 
of relatively large Einstein A. In third position, and with 
a sensitivity ratio of about 3, lies CH3OH, which also 
presents a depletion radius slightly larger than average, 
sizable A coefficients, and intensities similar to those of 
the SO lines. Interestingly, the molecule with largest de- 
pletion radius, C 18 0, has a sensitivity factor marginally 
better than 2, which results from a combination of a low 
dipole moment and moderate optical depth. The thinner 
C 17 (not shown) is significantly more sensitive, although 
still suffers from contamination by the extended cloud. SO, 
C2S, and CH3OH seem therefore the top three choices for 
any systematic search of young, undepleted cores. 



5.6. Origin of the line emission. Contribution functions 

Our L1498 and L1517B models reproduce the observed 
molecular emission, so they are expected to provide good 
approximations to the core internal excitation and radi- 
ation transfer. We can therefore use them to investigate 
how the line emission of different molecules is produced 
inside the core, and how the emission propagates and is 
distorted by optical depth effects as it travels towards the 
outside. Understanding these processes is critical when 
interpreting observations of similar systems for which a 
full radiative transfer analysis is not available, and can 
only be done by examining the internal properties of well- 
calibrated models. 

The general problem of line formation in a core is 
complex because of the non linear nature of the radiative 
transfer equation. In this section, we concentrate on the 
question of how the different layers of a core contribute 
to the emerging intensity, and on how different molecules 
can sample (or miss) the internal structure of the core. 
To quantify the discussion, we make use of the "contribu- 
tion function" (C F) commonly used in the stu dy of stellar 
atmospheres (e.g., iGrav hooA iMagain Ill98fih . This func- 
tion is derived from the formal solution of the radiative 
transfer equation. According to this solution, the emerg- 
ing intensity from the core at a given angle and frequency 
(without the background contribution) is 

Il ore = r S{r)e- T dr= f m j(r)e- T dl, 
Jo Jo 

where S is the source function, r the optical depth mea- 
sured from the core surface inwards, r m is the total core 
optical depth, j is the emissivity, dl is the line-of-sight el- 
ement, and l m is the full line-of-sight length of the core. 
In this notation, the contribution function in space units 
is 

V(!)=3(r) e- T[X) 

(note that stellar atmospheres texts commonly define the 
CF as a function of optical depth) . 

The emergent core intensity is the line-of-sight inte- 
gral of so the CF is true to its name in the sense 
that it measures the contribution of a given core element 
to the observed line intensity. Strictly speaking, the CF 
depends on wavelength because different regions of the 
core may contribute to different parts of the line profile. 
However, because of the small velocity gradients in L1498 
and L1517B (paper 1), and for the sake of simplicity, here 
we will only deal with the frequency-integrated CF, which 
measures the contribution of a line-of-sight element to the 
integrated line profile. Also for simplicity, we will only 
study the central line of sight of the core, although the 
generalization of the CF for non-zero impact parameters 
is straightforward. With these assumptions, it can be eas- 
ily shown that the CF of an optically thin line (where 
c & = 3) i n an isothermal core has two simple limits. If 
LTE applies ("high density"), the CF is proportional to 
the density n(r), so the emergent intensity is proportional 
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Fig. 10. Normalized contribution functions of representative species in L1517B (the observer is at large radius, i.e., 
to the far right of each panel). The red solid and dotted lines are the contribution functions, while the black dashed 
lines indicate the normalized density profile shown for reference (the normalized square of the density profile is also 
shown for CH3OH). In the best fit models (right column), the solid line corresponds to our simple abundance step 
function, while the d otted lines are the contribution functions expected for the abundance profiles predicted by the 
lAikawa et al. I l|2005l) models (scaled appropriately to fit the L1517B data). In each panel, the vertical dashed lines 
indicates the radius at which half of the total emission is reached. Note how in the moderately thick C 18 0(2-l) and 
extremely thick HCO + (1-0) lines most of the emission arises from the front outer layer of the core. 



to the column density. If A u ; 3> C u i (where A u i is the 
Einstein A coefficient and C u i is the collision coefficient, 
"low density" case), the CF is proportional to n(r) 2 , as 
each emerging photon is the result of a collision (which 
has a probability proportional to n(r) 2 ). In this case, the 
emergent intensity is proportional to the neutral analog 
of the "emission measure." Monte Carlo tests using differ- 
ent species and transitions under realistic core conditions 
show that the CF for a thin line usually lies somewhere 
between the above two limits. 

To illustrate the variety of CFs found during our radia- 
tive transfer modeling, we present in Figure 10 a series of 
normalized CFs for different species in L1517B (plots for 
L1498 are similar). Each panel represents a cut along the 
central impact parameter of the core, has the density peak 
at r — 0, and assumes that the observer is located to the 
far right of the plot (at r — > 00); a normalized density pro- 
file in dashed lines indicates the LTE optically thin limit 
of the CF. Constant abundance models are presented in 



the left column and best-fit abundance models appear on 
the right; we first discuss the constant abundance case. 

As Figure 10 shows, the constant abundance model 
for N 2 H+(l-0) (and that of NH 3 , not presented) has a CF 
that closely follows the density profile, with a a slight shift 
to positive radius due to optical depth effects (see below). 
This good behavior of the CF shows that the N 2 H+(l-0) 
emission responds linearly to density and therefore traces 
faithfully the core structure. To quantify this property, we 
compare the fraction of integrated CF that arises from the 
"inner core" (as defined within the half maximum density 
radius) to the fraction of gas column density contained in 
the same region (0.62 for L1517B). For both the N2H" 1 " 
and NH3 constant abundance models, we find that the 
inner core contributes to the emission a percentage indis- 
tinguishable from the ideal 62% value. In contrast, the CF 
of CH 3 OH(3o-2 ) in a constant abundance model (Fig. 10, 
second left panel) follows closely the n(r) 2 curve, and is 
therefore significantly biassed toward the high density gas. 
For this line, the inner core contributes 85% to the emer- 
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gent emission, which is almost 40% more than the column 
density fraction. 

The CF for the C 18 constant abundance model shows 
signs of non negligible optical depth in the form of a for- 
ward shift of the peak. This shift indicates that most of the 
observed C 18 0(2-l) emission arises from the front part of 
the core, and that the emission from the back is heav- 
ily attenuated as it travels through the core. In fact, only 
15% of the emerging radiation originates from the backside 
(r < 0), and the inner core contributes with only 38% of 
the emission (a 40% decrease from the ideal value). C 18 0, 
therefore, is a biassed tracer of the core emission despite 
its LTE excitation, although most of these problems can 
be minimized using C 17 transitions (but see below for 
models with depletion). As a final example of a constant 
abundance model, we consider the heavily self absorbed 
HCO + (1-0) line, although similar results are obtained for 
CS(2-1). The emission from this line is so dominated by 
the outer core that half of it originates from radii larger 
than 90", and the inner core only contributes 17% of the 
emergent intensity (23% in the case of CS(2-1)). These 
numbers illustrate the difficulty of modelling self absorbed 
lines like HCO+(1-0), CS(2-1), or HCN(l-O): the region 
that dominates the line emission is poorly constrained by 
the continuum observations that are the basis of our den- 
sity profile, and its structure is simply guessed using an 
extrapolation of the central core parameters. These very 
thick lines are therefore excellent tracers of the outer ve- 
locity fields (like infall), but they become poor tools for 
deriving any line of sight structure, even of the infall ve- 
locity, unless there is a very accurate description of the 
outer core and the surrounding cloud. 

The presence of depletion introduces a substantial dis- 
tortion in the CF. This is illustrated by the right-hand 
panels of Figure 10, which are based on the best-fit abun- 
dance profiles of section 4. The N2H + CF is of course the 
same as in the constant abundance case because no abun- 
dance gradient was found for this molecule. The CH3OH, 
C 18 0, and HCO + CFs, on the other hand, present sharp 
drops near the core center that correspond to the abun- 
dance holes found for these species. Because of these drops, 
the CF in the inner core vanishes or becomes so small that 
a negligible part of the emergent emission arises from the 
dense gas: 5% in HCO + (1-0) and less than 1% in both 
CH 3 OH(3o-2 ) and C 18 0(2-l). The spectra towards the 
core center for these species is therefore overwhelmed by 
gas outside the core, mostly from an innermost ring if the 
emission is not very thick (the case of CH3OH and C 18 0) 
or from the front part of the core in the thicker HCO + (l- 
0) line. Although it can be argued that the abundance 
drops in our models are artificially steep because of the 
use of a step function, theoretical considerations of core 
chemistry still predict very sharp abundance decreases. To 
test the effect of smoother, but realistic abundance gra- 
dients, we have produced a se t of CFs using the a = 1.1 
model of lAikawa et al. I l|2005l) corrected by the appropri- 
ate scaling factors to match the observed radial profile. 
The new CFs (shown as dotted lines in the right hand 



panels of Fig. 10) again present very steep holes toward 
the core centers and predict inner-core contributions to 
the emerging intensity of 15% (CH3OH) and less (C 18 
and HCO + ). The inner core contribution to the emission 
for these and other depleted species is therefore a minor 
fraction of the total emergent intensity. 

The above small sample of CFs illustrates some of 
the basic phenomenology of line formation in cores like 
L1517B. As the plots show, even in the simple case of 
constant abundance, the emission is not generated at a 
"t = 1 surface," but arises from a broad range of radii 
that in the simplest case (thin LTE) has a weight that 
mimics the core density distribution, and that in a thick 
line is strongly biassed toward the front part of the core. 
When realistic abundance profiles are considered, if the 
species is sensitive to depletion, all hope of tracing the in- 
ner core is lost no matter how much spatial resolution is 
achieved. Only depletion-resistant species like N2H" 1 " and 
NH3 trace the inner gas, at least to densities of a few 
10 5 cm~ 3 . Among these tracers, N2H + (l-0) seems spe- 
cially reliable because of a combination of constant abun- 
dance, reasonably low optical depth (helped by the hyper- 
fine structure), and well-behaved excitation. Ignoring the 
NLTE ratios between the J=l— hyperfine components (a 
10% effect most likely due to a population redistribution), 
this line consistently produces a CF very close in shape 
to the density profile. This is somewhat surprising given 
the significant radial gradient in the excitation tempera- 
ture of this line (more than a factor of 2 change between 
center and edge in L1517B), and it results from a popula- 
tion of the J=l level that follows closely the density law 
(the population of higher N2H + levels falls faster than the 
density law, and the CFs of higher transitions is closer 
to n(r) 2 ). This fortunate behavior of the N2H" 1 " molecule 
explains the striking similarity between the mm contin- 
uum and N2H + (l-0) maps in many st arless cores (e.g, 
Caselli et aflll999MTafalla et al. I l2002t iBacmann et al. I 
20021 ICransi et a,l.ll2005[ " this work), and makes the 1-0 
line the ideal choice for studying core interiors as long as 
its optical depth remains low. 

6. Conclusions 

We have presented observations of 13 molecular species 
(plus a number of isotopologues) toward the Taurus- 
Auriga starless cores L1498 and L1517B. Combining 
Monte Carlo radiative transfer modeling with a physical 
description of the cores, we have derived a self-consistent 
set of molecular abundance profiles. From the analysis of 
these profiles, we have reached the following main conclu- 
sions. 

1. Most abundance profiles can be described with sim- 
ple step functions having a constant value in the outer 
layers and a central hole. In both L1498 and L1517B, 
we find central abundance drops for CO, CS, CH 3 OH, 
SO, C3H2, HC 3 N, C 2 S, HCN, H 2 CO, HCO+, and DCO+, 
while N2H" 1 " and NH3 are found to remain in the gas phase 
even at the highest densities. The central abundance drops 
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suggest that most species freeze out onto the cold dust 
grains at the centers of the two cores. 

2. For most species, the abundances estimated in the 
outer layers of L1498 and L1517B agree with each other 
within a factor of 2, while the central abundance holes are 
systematically larger in L1498 by a factor of about 1.5. 
This suggests that the two cores have contracted from 
gas of similar composition, but that they are following 
different contraction histories or that they are at different 
evolutionary stages. 

3. Not all species with depletion present central holes 
of the same size. Molecules like DCO+, HCN, and HC 3 N 
systematically have smaller than average central holes, in- 
dicating that they survive in the gas phase up to slightly 
higher densities than the rest of the depleting molecules. 
Species like SO, C2S, CH3OH, and probably CO, on the 
other hand, present larger than average abundance holes, 
and seem therefore to deplete at lower densities. 

4. LI 498 and L1517B seem to have chemical composi- 
tions similar to other Taurus cores like TMC1, although 
no other multi-molecular survey of core chemistry that 
includes depletion is available for a detailed comparison. 
From a review of different abundance determinations in 
the literature, it seems that the chemical composition of 
L1498 and L1517B is rather typical of the Taurus starless 
core population. The abundances determined in this work 
can therefore be used as a reference for future modeling. 

5. A comparison with th e recent (and most up to date) 
model of core chemistry bv lAikawa et al. 1 J2005) shows a 
combination of successes and failures. For most molecules, 
the model predictions need to be corrected by factors of 
5 or even 10 in order to match the observed intensities. 
The model, on the other hand, predicts reasonably well 
(factor of 2) the size of the central depletion hole for most 
molecules. However, it fails predicting the observed variety 
of hole sizes among molecules. 

6. By comparing observed central intensities and inten- 
sities predicted by constant abundance models, we have 
studied the relative sensitivity of the different molecules 
and transitions to the presence of central depletion. We 
find that SO and C2S are the most sensitive indicators 
of molecular depletion, and therefore constitute the best 
choices when searching for very young cores. 

7. We have used radiative transfer modeling and re- 
alistic abundance profiles to study how molecular lines 
originate at the core interior. We find that the "contribu- 
tion function" provides an excellent tool to visualize and 
quantify the formation of the line profile from the differ- 
ent parts of a core. Using this function, we show that the 
N2H + (l-0) line is the most faithful tracer of the internal 
core structure as long as it remains optically thin. Species 
sensitive to freeze out only trace the edges of the central 
depletion hole (if thin) or the outer front part of the core 
(if thick). 
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Appendix A: Molecular parameters and 
convergence tests 

C2S. Dicarbon monosulfide is a linear radical with two un- 
paired electrons (electron spin S=l), and therefore fine struc- 
ture in its spectrum. Its rotational levels are characterized 
by N and J, the quantum numbers of rotational angular 
momentum and total a ngular momentum, respectively (e.g., 
IWolkovitch et al. Ill997l) . We model its radiative transfer us- 
ing the C2S molec ular parameters from the JPL catalog 
dPickett et al. Il l998l) to g ether with the C2S-H2 collision rates 
from IWolltcvi^ch"etal~l Jl997f) (downward rates for 10K to- 
gether with upward rates derived from detailed balance). We 
include all energy states up to approximately 60 K, which re- 
sults in a total of 37 levels and 79 transitions. 

CH3OH. Methanol is a slig htly asymmetric top with a hin- 
dered internal rotation llTownes fc Schawlow Il955l). It consists 
of two symmetry species denoted A and E ("e.g.. iLees II1973T) 
whose rotation energy levels are characterized by their total 
angular momentum J and its component along the near sym- 
metry axis k (an additional + or — is required for the A species 
if k^ 0). For our calculations, we have used t he energy and 
line s trength determinations in the JPL catalog iPickett et al. I 
1998 : ) together with the recent CILOH-para H2 collision rates 
of IPottage et af~l J2004T) (downward rates for 10K together 
with upward rates derived from detailed balance). For both 
species, we include all energy states up to approximately 70 
K, which results in a total of 31 levels and 92 transitions for 
E-CH3OH, and 27 levels and 66 transitions for A-CH3OH. 

C-C3H2. Cyclopropenylidene, a ring molecul e, is an oblate 
asym metric top with ortho and para species iVrtilek et al. I 
1987). Its rotational levels are characterized by the total angu- 
lar momentum J and by K_i and K+i, its two projections 
on the symme try axis for the prolate an d oblate symmet- 
ric top limits llTownes fc Schawlow I Il955l) . For our calcula- 
tions, we have used the C3H2 parameters in the JPL catalog 
dPickett et al. I l998l) tog ether w i th the C3H2-He collision rates 
for 10 K from Ia very fc Green I dl989l) (multiplied by 2 1/2 to 
simulate collisions with H2, as recommended by the authors). 
Only transitions of ortho C3H2 were observed, so our modeling 
is restricted to that species. This modeling includes all states 
with available collision rates (up to E u = 41 K): 16 levels and 
32 transitions. 

H2CO. Formaldehyde is a slightly asymmetric prolate ro- 
tor with ortho and para species, and its rotation levels are 
characterized (like those of C3H2) by J, K_i, and K+i. For 
our calculations, we have used the H2CO parameters from the 
on-line Cologne Da tabase for Molecular Spectroscopy (CDMS, 
iMiiller et al. 1200 ill, t ogether with the 10 K H 2 CO-He collision 
rates of |GjeeriT<ll99 ll) . These collision rates were multiplied by 
2.2 to convert them into H2CO-H2 rates, as recommended by 
iGreen I <ll99ll) . We only observed lines of ortho H2CO, so only 
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this species has been modeled. For this, we have taken into 
account all energy states with energy below 115 K, resulting 
in a calculation with 14 levels and 19 transitions. 

HC3N. Cyanoacetylene is a linear molecule. Its rotational 
levels have hyperfine (hf) structure due to the nuclear spin 
of the nitrogen atom, so they are characterized by the quan- 
tum numbers J and F. For our radiative t ransfer calculation 
we hav e used the molecular parameters fromjLaffertvj^^ova^ 
W7& ) and the collision rates with He of iGreen fc Chapman I 
1978) (multiplied by 2 1//2 to simulate collisions with H2, as 
recommended by the authors). The calculation treats the hf 
structure in a simplified manner by solving the equations of 
statistical equilibrium for the pure rotational levels after di- 
viding by 3 the effect of trapping, to simulate the decrease of 
optical depth due to the three-fold hf splitting (we find this to 
be a minor effect) . The resulting population of each rotational 
level is then divided between the hf su b levels assuming LTE 
to compute the emergent spectrum fsee lTafalla et al~ll2002l for 
a discussion of a similar approach with N2H + ). All 24 levels 
up to an energy of 130 K are considered (23 transitions). 

HCO + and isotopologues. The formyl cation HCO + 
and its isotopologues H 13 CO+ and DCO+ are linear ions. We 
have modeled their radiative transfer using the parameters in 
the JPL catalog togeth er with the 10 K HCO + -H2 collision 
rates of iFlowerl 1^999). Our calculation considers a set of 7 
levels (6 transitions), up to an equivalent energy of 120 K. To 
take into account the small additional broadening due to hy- 
perfine structure in H 13 CO + , we have artificially broadened 
this l ine by an extra 0.133 km s _1 (see ISchmid-Burgk et al.H 
l2004ft . 

HCN and H U CN. Hydrogen cyanide is a linear molecule 
with hyperfine structure due to the nuclear spin of the nitro- 
gen atom. We have modeled its radiative transfer using the 
molecular parameter s in the CDMS together w ith the HCN- 
He collision rates of iMonteiro fc St utzki I i l986T) ( f or J u p to 
4) and the HCN- He rates of IGreen fc Thaddeusl < 197-1) (re- 
maining levels up to J= 7). A factor of 2 was used to convert 
the He collision rates into H2 rates. The hyperfine structure of 
this species (and that of its isotopologue) has been considered 
explicitly in the calculations, although overlap effects were ig- 
nored and the collision r ates between hf sub l evels with J> 4 
were estimated from the lGreen fc Th ackdeusJ 1^7^) rates as- 
sumi ng independence of the F number IIGuilloteau fc Baudrv I 
198 ll) . A total of 22 levels and 39 transitions were included, 
corresponding to all levels with energy lower than 119 K (up 
to J=7). 

SO. Sulfur monoxide is a diatomic molecule with two un- 
paired electrons and therefore a fine structure similar to that of 
C2S. We have modeled its radiative transfer using the molecu- 
lar param eters in the CDMS t ogether with the SO-He collision 
rates from lLiaue et al~l i2005T) . Following the authors's recom- 
mendations, a factor of 1.38 was used to convert the SO-He 
collision rates into SO-H2 rates. All 24 levels up to an energy 
of 85 K and their corresponding 72 transitions were included 
in the model. 

Convergence tests. When solving the radiative transfer, 
we divide the core in 200 equally spaced shells. We normally 
use 2000 photon packages to simulate the radiation, and allow 
40 iterations for the so lution to conv erge (no radiation refer- 
ence field was used, see lBernes 1119791 for details on the code). 
To check the convergence of the solution, several tests were 
carried out. Independence on the initial level population was 
tested by comparing the results from runs that started with 




Fig. B.l. Spectra towards L1498 illustrating the presence 
of a red gas component. The bottom (H 13 CN) line shows 
that the LSR velocity of the cloud is about 7.8 km s _1 
(vertical line). The HCN line (next up) shows a deep red 
self- absorption and a low level red wing, while the CO 
isotopologue lines show an increasing contribution of the 
red component with abundance. The coincidence of the 
red component with the absorption in HCN (also HCO + 
and CS) suggests that the component lies in front of the 
L1498 core. 

Boltzmann equilibrium at 10 K (kinetic temperature) and at 
2.7 K (cosmic background); no differences in the final solution 
were found. Independence on the internal Monte Carlo param- 
eters was tested by repeating the best fit solutions using 5000 
photon packages and 300 iterations. Again, no appreciable dif- 
ferences were found in the predicted intensities. 

Appendix B: An additional gas component toward 
the L1498 core 

The line of sight toward the L1498 core is complex. In addi- 
tion to the narrow-line component that corresponds to the core 
studied in this paper, spectra of abundant species like 12 CO 
and 13 CO show extra emission toward the red extending up to 
about Vlsr = 9.5 km s _1 , which is 1.7 km s _1 higher than 
the core velocity (Fig. B.l). This add itional component h as a 
complex spatial distribution (see also lKuiper et ai~lll996T) . At 
low velocities (~ 8.2 km s - 1 ), it lies toward the south west of 
the core (as also found bv lLemme et al."lll995l while at high 
velocities (~ 9 km s _1 ), it appears toward the NW. This distri- 
bution around the L1498 core suggests some relation with the 
dense gas, although the velocity difference implies that the red 
gas is not bound to the core. From the C 18 O(l-0) spectrum 
in Fig. B.l and assuming an excitation temperature of 5-10 
K together with a "s tandard" C ls O abundance of 1.7 10 -7 
lFrerkin~ et al. I (I1982D . we deduce an approximate H2 column 
density for this component of 2 10 21 cm -2 (less than 10% of 
the core column density). 
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Although the red component does not affect the narrow 
optically thin lines used in our abundance determinations, it 
can distort the shape of very thick spectra like HCN(l-O), 
shown in Fig. B.l. These spectra were also modelled in section 
4 by assuming standard isotopic ratios, and they sometimes re- 
quired anomalous abundances in the outer core (section 4.2). 
As stressed before, these optically thick lines were only fitted 
to prove self consistency, and not to determine the conditions 
of the outer core layers. To model these lines, however, some 
assumptions had to be made about the red component, and 
a main one was its location along the line of sight with re- 
spect to the L1498 core. A foreground position was preferred 
based on hints from several optically thick spectra, like the 
HCN(JF=12-01) line shown in Fig. B.l. This spectrum shows 
a deep self absorption that removes completely the red part of 
the line and replaces it by a low-level plateau that continues 
smoothly into the red-component range (from Vlsr 8 to 9 km 
s _1 ). Such a feature is best understood if the red component 
contributes to the absorption of the core emission and there- 
fore arises from gas located in front of the core. To model such 
a component, we have red-shifted the front outer core layers of 
L1498 ("envelope" in paper 1) by 0.35 km s _1 and broadened 
them to AV = 0.625 km s -1 (note that because of a typo, 
Table 3 in paper 1 gives the wrong envelope velocity, which 
should be 8.15 km s _1 ); such a simple parameterization repro- 
duces approximately the shape of most optically thick lines. 

Appendix C: Comparison with published data 

We have compared our L1498 model predictions with data 
for this core availab le in the literature (no s imilar data were 
found for L1517B). IWolkovitch et al~l (119971) have presented 
observations of L1498 in three C2S transitions, JN=21-01 (22 
GHz), JN=43-32 (45 GHz), and JN=87-76 (94 GHz), made 
with the telescopes of NASA's Deep Space Network and the 
FCRAO, and having spatial resolutions of about 50". From 
the linewidths and peak intensities reported by these authors 
in their Table 2, we have reconstructed the line integrated in- 
tensities for all positions available, and with them, we have gen- 
erated the radial profile of emission for each transition. When 
we compare these profiles with the prediction from the best 
abundance model derived in section 4.2, we find the excellent 
match illustrated by the three top panels in Fig. C.l. This 
match shows that our model predictions are ful ly consistent 
with the observations of IWolkovitch et al. I il997t) . 

Other L1498 da ta of interest for our study are those of 
lYoung et aTl <2004h . who have observed L1498 in H 2 CO(ln- 
lio) (6 cm) with the Arecibo telescope and in H2CO(3i2-2n) 
(1.3 mm) with the CSO (spatial resolutions of 60" and 32"). 
To reconstruct the integrated intensities of these lines, we have 
used the values reported by the authors in their Table 2 (1.3 
mm line) and data kindly provided by Kaisa Young and Jeong- 
Eun Lee (6 cm line). With them, we have generated a set of 
radial profiles, and in Figure C.l we compare them with the 
predictions from the best fit model of section 4.2. As the figure 
shows, the model reproduces well the 3i2-2n emission, which 
is sensitive to the high density gas inside the core. The model, 
however, fails to reproduce the observed intensity of the I11- 
lio transition, although it agrees with the data in producing 
an absorption line against the cosmic background. The poor fit 
occurs because the absorption depends sensitively on the low 
density gas outside the core, and this gas is not well parame- 
terized by our extrapolation of the inner core model (section 
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Fig. C.l. Comparison of L1498 data taken from the lit- 
erature (solid squares) with the predictions of the best 
fit model of section 4.2 (red lines). The top three pan- 
els show that the model fits well the C?S radial pro - 
files generated from the data of lWolkovitch et al. 1 (ll997lL 
The bottom two pane ls show that the mo del fits well the 
H 2 CO(3i2-2n) line of lYoung et aTl ll2004h . but fails to fit 
the H2CO(ln~lio) line (from the same authors), which 
appears in absorption against the cosmic background. To 
illustrate the sensitivity of the ln-lio absorption to the 
core low density g model with a backside envelope is 
shown in blue dashed lines (see text). 

4.2). To illustrate this sensitivity to the low density gas, we 
have added to the L1498 model an envelope in the back that 
has the same properties as the front envelope used to simulate 
the red shifted component (appendix B). This back envelope 
doubles the ln-lio absorption, and brings the model predic- 
tion closer to the observations (blue dashed lines in Fig. C.l). 
Such strong sensitivity to the low density gas shows that the 
central core is close to invisible in the ln-lio transition (note 
the flat radial profile), and that improving the fit requires a 
re-parameterization of the envelope, not the core. Given the 
few observational constrains of the low density gas, such a pa- 
rameterization lies outside the scope of this work. 
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